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We develop a forward algorithm for solving and finding planning horizons for the infinite horizon version of the deterministic 
production smoothing problem without inventory. This model approximates the real world situation for highly obsolescent or 
perishable commodities such as newspapers and fresh produce. Computational results show that the algorithm is linear in problem 
length while linear programming is at least quadratic. 


1. Introduction 


The finite horizon production smoothing problem with linear costs has been considered by a number of 
authors [1,5,6,8,9~11,14]. Here we consider the special case in which neither stockouts nor inventories are 
allowed. This model approximates the real world situation for highly obsolescent or perishable commodi- 
ties such as newspapers and fresh produce. 

We consider the infinite horizon version of the problem which, after Lundin and Morton [12] and 
Morton [13], involves solving longer and longer finite horizon problems with a ‘forward algorithm’. Only 
the optimal decisions for a first block of periods termed the ‘planning horizon’ are desired, since the model 
is actually re-solved periodically to incorporate improved forecast information (a rolling horizon proce- 
dute). As soon as the decisions in the planning horizon can be guaranteed not to change for any longer 
finite horizon problem a ‘forecast horizon’ is said to be achieved, and computation stops. If for some 
reason a forecast horizon is not achieved, computation stops after solving some ‘sufficiently long’ finite 
horizon problem. 

The formal model is presented in Section 2. Although for a fixed horizon length it may be solved by 
linear programming, it is not efficient to construct a forward algorithm simply by solving a sequence of 
LP’s. The special structure of the problem is exploited to develop an efficient forward algorithm for the 
problem in Section 3. In Section 4, planning horizon procedures are developed to provide stopping rules for 
the algorithm. A twenty period example is solved in Section 5. Finally in Section 6, experimental results are 
presented which demonstrate that the algorithm is superior to linear programming, being about linear in 
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computational requirements in problem length. Furthermore planning horizons are in practice found very 
quickly. Proofs for the theorems and lemmas appear in the Appendix. 


2. The model 


Let an infinite sequence of demands D,, D;, ... be specified. The workforce in period ¢ is W,. overtime 
is V.. undertime is U.. The number of workers hired at the beginning of period fis H,, the number fired is 
F.. The nonnegative costs of overtime, undertime, hiring, and firing per unit are v, u, A, and f respectively. 
The T period problem (for T= 0, 1, 2, ...) is given by 


T 
(P,) minimize YY (uUt+eV.+AH,+fE)+hH,,,4+fF,.,. (1) 
f= 
subject to 
WW. ,-H,+F,=0, t=1..... T+, (2) 
W+V-U=D, P= lots T. (3) 
W,, 1 = W.. (4) 
W.H,.F.U.V,20 where W,, W, are given constants, (5) 
D,, W,. W, > 0 integer. (6) 
Since demand is to be met for all ¢= 1,.... T we do not include the production costs in the objective 


function (1). 
Clearly (P,-) is a linear program, and may be solved by standard techniques such as the simplex method. 
We adopt linear programming terminology for the following results. 


Lemma 1. For any T= 0. 1..... there exists a basic feasible integer solution to (P;). 


Lemma 2. /f a basic feasible solution to (P,-) satisfies 


W,> mini H.W. min {D);. =) eer T+. (7) 


then every basic feasible solution found by the simplex method will satisfy (7). 


Note that by Lemma 2 we could drop the nonnegativity requirements on W,, since if we start with a 
solution satisfying (7), we can never pivot to one which violates (7). 

By inspection of the matrix, A, which represents the constraint set of problem (P,-), we see that the 
complementarity conditions 


H.-F =0, (8) 
V-U=0 (9) 


are satisfied by any basic feasible solution, otherwise the basis will contain linearly dependent columns. 
Define problem (P;’ ) to be identical to (P,.). except that all variables are primed and 


Wi= Wt, Wl =W.+N, Di/=D,+N, t=1...., T (10) 
where N is any nonnegative integer. Theorem | states that we can find an optimal solution of (P;) from 
one of (P;. ). 

Theorem 1. /f 
We de POs oe (11) 
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is an optimal solution to (Pf), then 
Wea We SNA PS PUT (12) 
is an optimal solution to (P,). 


So without loss of generality we may assume W,, W., D, > 0 for (P;), and that any initial basic feasible 
solution satisfies (7). 


Theorem 2. Jn every non-degenerate primal simplex pivot for (P+). exactly two hiring and / or firing variables 
change value. 


Theorem 2 hints at the procedure we discuss later in this paper. If we could avoid degenerate pivots by 
insuring that the indicated change occurs, we could save on needless basic updates. 


Theorem 3. Ail basic feasible solutions of (P;) are integer. 


So, we do not consider non-integer changes in H, and F, during a pivot. 
Since all costs are nonnegative, there exists an integer optimum to (P;). For any feasible solution S,. to 
(P;), and any pair of periods 1, j, 1 <i<j<T+ 1, define: 
k(i, 7, S;) = the number of periods V,> O,i<1t<j, 
<1<j, (13) 
n(i,j, S;) =the number of periods U. = V,=0,i <1 <j. 


m(i,j, S;) = the number of periods U, > 0, i 


Further define: 


X(i, j, S;) = vk — um — un, as the net benefit of adding one worker 

from period i toj — 1. (14) 
Y(i,j, S;) = um — vk — vn, as the net benefit of removing one worker 

from period i toj — 1. 


Whenever S; is deleted from the above definitions, we infer ‘the current solution’. 

We also define a number of pivor steps for improving solutions. The notation H, f F, 1 means: increase 
hiring in period / and firing in period j by an equal number of units until overtime in one period from i to 
jJ— 1s first driven to zero. 

The notations H,t H,|, F.lH,1, F,) #1 are defined similarly, except that the H/F increase or 
decrease is terminated also if an H/F variable reaches zero before an overtime variable does. Similarly, for 
driving undertime in one period from i to / — | to zero, we define the pivots: F, t Wt ETF Hl Fil, 
HV H,t. 


Table | 
Pivots 
Condition Pivot Unit cost reduction 

(i) H,>0,F 20, X(i,j)>(A+f) Ht Ft (A+ f)- X(i7) 
(ii) H,>0,H,>0, X(i,j)>0 H,1H,4 - X(i.j) 

(ili) F,>0,H,>0, X(i,j)> -(h+f) Fila, —(ht+ f)- Xi) 
(iv) F>0,F, 20, X(i,7)>0 FLFt — X(i.j) 

(vy) F20,H, 20, ¥Q.j)>(h+f) Fth,t (A+ f)-YC.Y) 
(vi) F,20,F >0, ¥Cs)>0 FTEs -YQ.7) 

(vii) H,>0.F>0, Yj) > —(h+f) HVE -—(h4t f)- YC.) 


(viii) H,>0,H, 20, YUi,j)>0 H,.H,?T -—Y(4./) 
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Lemma 3. In each of the eight cases given in Table \, for a current policy S, and periods i and j, the indicated 
pivot will lead to a strictly improved solution 


It can be easily shown that the other eight pivots possible can only lead to an infeasible solution. So, we 
have narrowed our scope to only eight possible improving pivots. 


The following result allows us to determine when a solution to (P,-) is optimal. 


Theorem 4. S,. is an optimal solution if and only if for no pair of periods \ <i <j < T+ 1 do any of the eight 
conditions of Table 1 hold. 


This is, in effect, a pricing check on non-basic variables. If any of conditions (i)-(vili) hold, then an 
improving pivot step can be performed. We next state two conditions that occur in optimal solutions. 


Lemma 4, For an optimal plan S¥, D, > W, implies F. = 0: D, < W, implies H, = 0. 
That is, don’t fire or hire needlessly. 
Define 


r(S,) = max{i| X(i, T+ 1) > X(i’. T+ 1) for all i’}, 
> 


1S 
s(S;) = max{i] ¥(i, T+ 1) > YC“, T+ 1) for all i’). dy) 


In other words, the pair (7, T + 1) is chosen to have a larger potential for decreasing overtime than any 
other pair (i, T+ 1); the pair (s, T+ 1) is chosen to have a larger potential for decreasing undertime than 
any other pair (/, T+ 1). 

In the following r = r(S#) and s = s(S#) 


Lemma 5. For an optimal plan S¥, X(r. T + 1)>0 implies F,=0: Y(s. T+ 1)> 0 implies H, = 0. 


Again, don’t fire or hire needlessly. 
Further define 


Pp, =max(t|H,>0,3t’>1: F.>0} (0, if undefined), 


16 

q, = max(t| F,>0,30 >t: H,>0} (0. if undefined). uly 
p =min{t|H,>0,¢>q,) (0 if undefined). 

=min({t|F,>0.t>p,) (0 if undefined) (17) 


w = max{ p.q). 


In other words p is the first hiring period following the last firing period q,: and q is the first firing 
period following the last hiring period, p,. 
We also define sets 


P=({t|H,>0,1>p)}., Q=P(t|F,>0.1>q)}. (18) 


The set P contains periods for which hiring occurs following the last firing period, similarly for set Q. Note 
that either or both of these sets may be empty. In the next section, we describe a Forward Algorithm which 
uses the indices defined earlier and the sets P and @Q to solve (P;.). 


3. The Forward Algorithm 


The Forward Algorithm is outlined in Fig. | (see also Fig. 2). Each time the plan is changed, either by 
augmentation to a longer horizon plan. or pivoting to an improved plan. we verify that the period 
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1. Initialize by solving the 7 = 0 problem. 
2. Set T=7 +1. Given the optimal solution S#_, to the T—1 period problem, augment it to form a feasible solution S!. to the T 
period problem. 
3. Pivot as many times as is necessary to drive S!. to S¥. Three separate procedures are described depending on whether 
a. Dp = Wh, 
b. D; > W., 
ce. Dr < W,. 
4. Go to Step 2. 


Fig. 1. Forward Algorithm outline. 


indicators w = max{ p, q}, and r are (weakly) monotonically increasing. This provides strong computational 
savings in implementing the Forward Algorithm. It also provides the basis for the planning horizon 
procedure given in the next section. 

We turn next to the question of pivoting on the augmented plan. In Lemma 6 we show that H;,, or 
F,., must change value to improve the plan. 


Lemma 6. The only possible pivots on S} are for pairs of periods (i, T+ 1), forl <i< T+. 


Next, we restrict our attention to, at most, two possible pivoting conditions, depending on whether 
D,> W,, Dr < W. or Dp = W,. 


Lemma 7. The only possible (Lemma 3) pivoting conditions for S;.on (i, T + 1) are shown in Fig. 3. 
So, at least initially, only conditions (i) and (iv), or (v) and (viii) could occur. In the case D,.= W,, there 


are no pivot steps because the augmented solution is optimal. 
We next consider pivoting for the case D; > W,, that is in (S}.), there is overtime in the final period. 


1. Solution to the T =0 Problem ( Initialization) 


1 W=W. 

2.r=s=p=q=p,=49,=0. 

3. P=Q=®. 

4. a. If Wy > W,, then F, = W, — W., and by convention: g =1. Q = {1}, 


b. If Wo < W,, then H, = W. — Wo, and by convention: p = 1, P = (1). 
5. By convention, X(0, i) = 0 for all ¢ = 1, 2,... 


2. Augmentation of the T —\ Period Problem (S#_, to S;) 
. Set T=T +1. 


I 
2. All decision variables of the T — 1! period problem remain unchanged. 
3. Wr. = WW, (= W,). 
4. Hy, = Fy, = 0. 
5. 
o r(Sp_,) if X(r.T4+1)> X(T. T +1), 
r(Sp)= 
T otherwise, 
‘ . 
yee (CER if tees 1)>Y(7T.T+4+1). 
T otherwise. 


6. a. If Dp =W,, then Up = V, =0 
b. If Dp > W,, then V;= D, -—W,, U;=0 
c. If Dy < W,, then U; = W.- Dy, V7 =0 


Fig. 2. Details of Steps 1 and 2 of Fig. 1. 
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Case Possible pivoting condition 


a. D, =W, None S$}. is optimal 
b.D,>W. (id (iv) 
co D, <W, (v), (vill) 


Fig. 3. Possible pivoting conditions for S}. 


Pivoung S} to S$ when D, = W, 
|. All decision variables are unchanged. 
2. a. p. gq. Py. gq) are unchanged. 

b. sets P, Q are unchanged. 


Fig. 4. Elaboration of case (a), Fig. 3. 


Theorem 5. For pivoting conditions (1) and (iv) the following seven statements hold: 
(a) If for S} pivoting condition (iv) holds for (i, T + 1). then it holds for (q*. T + 1), where 


q* = max{i| X(i,. T+ 1) > XC. T+ 1). 1. EQ). (19) 


(b) If for S} pivoting condition (i) holds for (i, T + 1). then it holds for (r. T + 1). 
(c) If for S} only pivoting condition (iv) holds for (i, T + 1). then after executing pivot (iv) on (q*, T + 1), 
only pivoting condition (iv) may hold, and only for (i, T + 1). 
(d) If for S} only pivoting condition (i) holds for (i, T + 1), then after executing pivot (i) on (r, T + 1), only 
pivoting condition (1) may hold, and only for (i, T + 1). 
(e) If for S} both pivoting conditions (i) and (iv) hold for (i, T + 1) and (j, T + 1) then 
(1) if r<q* or ifr>q* and X(r, T+ 1)-(h+f)< X(q*, T+ 1), then after executing pivot (iv) on 
(q*. T + 1). only pivoting conditions (i) and (iv) may hold, and only for (i’, T + 1) and ( j’, T + 1)). 
(2) ifr >q* and X(r, T+ 1)-(h+f)> X(q*, T + 1), then after executing pivot (i) on (r, T + 1), only 
pivoting conditions (1) and (iv) may hold, and only for (i’, T + 1) and (j', T + 1). 
(f) Executing pivot (iv) by the above rules does not cause uny of w and r to decrease. 
(g) Executing pivot (i) by the above rules does not cause any of w and r decrease. 


The reason for the special pivoting rule for when r > g* is that a pivoting condition may be introduced if 
the wrong pivot is performed. Let’s consider performing the ‘wrong’ pivot. Suppose r > qg*. By definition 
q* © Q and hence g* > p. so H,=0. If X(r. T+ 1)-(h4+f)< X(q*, T+ 1), then since X(q*, T+ 1)= 
X(q*.r)+ X(r, T+ 1), we have X(q*.r)> —(h+f). So if we perform the type (i) pivot, the type (111) 
pivoting condition is introduced, since X(q*, r) is unaffected by pivots on (r, T + 1). Now if X(r, T+ 1)- 
(h+f)> X(q*, T+ 1)>0 and pivot (iv) is executed instead of (1). it is possible to generate a type (v) 
pivoting condition on (q*, 7). This condition does not occur when following the rules specified by the 
theorem. 

The theorem presented above essentially outlines the algorithm which we develop to solve such 
problems. Based on the theorem we present the following corollary without proof. 


Corollary 1. Starting from S$}. the augmented solution of S*. ,, executing the pivots as outlines in Theorem 
5(c)-(e) leads to S} in a finite number of steps. 


Although counterexamples can be found to show that s is not necessary even weakly monotonically 
increasing, it can easily be shown that: 


Lemma 8. /f D,.> W, and at least one type (i) or (iv) pivot is required to drive S; to S}. then no type (v) pivot 
will ever be required after k augmentations for (j, T+ 1 +k), where j<T+ 1. 
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1. Wf X(¢*, T+1)>0 and g* > r; or if X(q*, T+ 1)>0. g* <r. and X(r, 7 +1)-(h + f) < X(q*. T +1), perform a type (iv) pivot 
and update p, g, r, s, and sets P, Q. 
Repeat. 

2. Otherwise if x(g*, 7+1)>0 and q* <r, X(r,T+1)-(h +f) > X(q*. T +1), perform a type (i) pivot and update p. q, r,s, and 

sets P, Q. Return to Step. I. 

. if X(r,T +1) < (A +f), STOP, the current plan is optimal. 

. Otherwise, perform a type (i) pivot on (r, T + 1), and update p, q, r,s, and sets P, Q. 

5. Return to Step 3. 


aw 


Fig. 5. Pivoting S; to St when D, > W,. 


1. If Y(p*, 7 +1)>0 and p*>s; or if Y(p*, 7 +1)>0, p*<s, and ¥(s,7+1)-(A + f)< ¥(p*. T +1), perform a type (viii) pivot 
and update p, q, r, s, and sets P. Q. 
Repeat. 

2. Otherwise, if Y( p*, 7 +1) > Oand p*<s, Y(s,T+1)-(h+/f)2> Y(p*.T + 1), perform a type (v) pivot and update p. q, r,s. and 

sets P, Q. Return to Step 1. 

. If ¥(s, 7 +1) < (A +f), STOP. the current plan is optimal. 

. Otherwise, perform a type (v) pivot on (s, T +1), and update p. q, r. s, and sets P. Q. 

5. Return to Step 3. 


> wW 


Fig. 6. Pivoting S} to S# when D; < W,. 


So, even though s is not necessarily increasing with the performance of type (i) and (iv) pivots, we may 
arbitrarily set s = 7 + 1 after any type (i) or (iv) pivot is performed, since we know that no type (v) pivot is 
possible for s < T+ 1. 

We show the pivoting step of the algorithm for the case of overtime in the final period in Fig. 5. Note 
that since Step 3 may never be executed more than once for each value of r, termination is assured. 

We define p* analogously to q*. In Fig. 6, we give the pivoting procedure for the case of undertime in 
the final period without proof. 

The complexity of the Forward Algorithm is described below. 


Theorem 6. Through completion of the T period problem, the Forward Algorithm pivots no more than 
T(T — 1)/2 times; computation is no worse than O(T*). 


The algorithm is thus potentially much more computationally efficient than the conventional simplex 
method. Following each augmentation, there is a limited number of non-degenerate pivot steps to examine 
and execute. The pricing operation and ratio test are trivially executed. The ‘candidate list’ for changing 
hiring and/or firing variables is limited by the planning horizon. As we shall demonstrate in the 
computational results in Section 6, the Forward Algorithm is linear in problem length 7, while the simplex 
method is worse than quadratic (at least O(7?)). In addition, the Forward Algorithm format allows 
computation to stop as soon as a planning horizon is achieved. This is discussed in more detail in the next 
section. 


4. Planning horizon theorems 


We are now ready to develop planning horizon procedures. Recall the definitions of p, and g, from (16). 
They may alternately be defined as: 
p,(S;) = max{t| H,>0,1<q(S;)} (0 is undefined), 


20 
qi\(S;) = max({t]F>0,1<p(S;)} (0 is undefined), me 
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We further define 
zp=min{w.r(S*), 5(S*)}-1 (21) 
where by convention r (resp. s) is set to 7 + | after a type (v) or (vill) (resp. (i) or (iv)) pivot. 
We call periods | to z, the stable part of the T period problem since decisions in those periods are 
unaffected by increasing the problem length from T to T + 1. 


Lemma 9. /f z; > 1, then periods | to T form a forecast horizon for a planning horizon of periods \ to zy. 


So, the decisions made in periods | to z, are unaffected by the demand in period 7 + 1, T + 2. etc. 


Lemma 10. 
(a) r(S}) > 4,(S#). 
(b) s(S#)>p,(S}). 


That is, the index r exceeds the period in which firing occurred followed by hiring in S*. Similarly, s 
exceeds the last period in which hiring occurred followed by firing. We now state the Planning Horizon 
Theorem. 


Theorem 7. z; > max{ p,(S¥%), q,(S})) — 1: that is to sav, as soon as any T period optimal plan has periods i, j 
such that i <j and either 


(a) H,>0,F,>0 or 
(b)  F>0,H,>0, 


then a planning horizon has been obtained at least through period i — 1. 


This means that far from planning horizons being rare, not obtaining a planning horizon requires a 
rather special demand pattern for which repeated hiring and firings do not occur. 

For a fairly typical case in which hiring and firing occurs each year, a planning horizon will be obtained 
in between a horizon length of one to two years. This does not depend on any particular demand pattern, 
only on successive hirings and firings. 


5. An example 


In Tables 2 and 3, we have tabulated results of a 20 period problem. In Table 2, the first column, 1, is 
the ‘current’ problem length. The second column, D,, is the demand in period t. The next 6 columns are the 
indices following the pivot, or augmentation in the event that there is no pivot. The first index, r, is the best 
type (1) pivot, s, the best type (v) pivot, p is the furthest period in time to which we must look for a type 
(vill) pivot, g for a type (iv) pivot. Note that p = min{t|t© P}, q= min(t|t© Q}, and if P (resp. Q) is 
empty, then p (resp. q) is set to 0, denoted by a — in this table. Further we have w = max{ p, g} and the 
planning horizon z, = min{w. r,5)}— 1. The set P (resp. Q) is the set of candidates for a type (viii) (resp. 
(iv)) pivot. The pivot period is chosen by the rules detailed in Figs. 5 and 6. 

The planning horizon z,, indices p and q are undefined until period 3, when the first type (v) pivot, on 
(1.4), is performed. Then, set P is no longer empty and we have a valid planning horizon z, = | for the 
forecast horizon T = 3. When T = 6 it becomes ‘cost effective’ (¥(4,7) > 0) to perform a type (viii) pivot on 
(4.7), thus decreasing the workforce from period 4 to 6 in addition to the reduction from period 1 to 3. The 
same is true for the type (viii) pivot performed on (7,8) when T = 7. 

The type (i) pivot for when T = 9 forces s = 10 and allows no more type (v) or (viii) pivots on (7. T+ 1) 
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Table 2 
Tabular results of an example of the Forward Algorithm (u = 44, v = 50, h = 54, f = 43, Wy = 58. W. = 58) 
t D, Indices following augmentation and Sets Pivot Cost 
if necessary, pivot P Q type (is) before after 
r s P q w ns) pair Pivot pivot 
1 48 l 1 - - - - - - 440 
2 38 2 1 - - - - - - 1320 
3 23 4 2 4 - 4 1 {4} - (v) (1.4) 2860 2510 
4 97 4 2 4 - 4 1 {4} - 4460 
5 30 4 2 4 - 4 1 {4) - 5692 
6 49 7 2 4 - 4 ] {4,7} = (vill) (4.7) 6088 5746 
7 20 8 2 4 - 4 1 {4,8} - (viii) (7.8) 7418 7022 
8 59 8 2 4 - 4 l (4,8) - 7072 
9 98 9 10 - 10 10 8 - {10} (i) (8,10) 9072 9069 
10 $7 9 10 - 10 10 8 . {10} 9113 
1! 33 9 10 - 10 10 8 - {10} 10213 
12 61 12 10 - 10 10 9 - (10) 10363 
13 76 12 14 - 14 14 NH - {14) (iv) (10.14) 11263 11251 
13 14 - 14 14 12 - {14} (i) (12,14) 11251 11245 
14 60 13 15 - 14 14 12 ” {14,15} (iv) (14,15) 11345 11245 
15 59 13 16 - 14 14 12 - (14,15,16} (iv) (15,16) 11295 11245 
16 62 13 17 = 14 14 12 - (14,15,17} (iv) (16.17) 11445 11395 
16 17 - 14 14 13 - {14,17} (iv) (15.17) 11395 11389 
17 18 16 17 - 14 14 13 - (14,17} 13149 
18 15 16 17 - 14 14 13 - {14,17} 15041 
19 15 20 18 20 - 20 17 {20} - (v) (17,20) 16933 15533 
20 24 21 18 20 ~ 20 17 {20,21} - (viii) (20,21) 17029 15533 
Table 3 
Optimal solution of the example 
‘ D, W, v, U, H, F 
1 48 48 10 
2 38 48 10 
3 23 48 25 
4 97 49 48 1 
5 30 49 19 
6 49 49 
7 20 49 29 
8 59 59 10 
9 98 59 39 
10 57 59 2 
al 33 59 26 
12 6] 61 2 
13 76 61 15 
14 60 60 J 
15 59 60 1 
16 62 60 2 
17 18 18 42 
18 15 18 3 
19 15 18 3 
20 24 24 6 
21 - 58 - - 34 
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where j < 10 (see Lemma 8). Also. since r = 9 and q = 10, the horizon jumps to zy = 8. 

In Table 3, we show the optimal solution for the 20 period forecast horizon. Note that zeros are 
suppressed and a—represents an undefined quantity. Due to Theorem 7 we have z,, = 17 > max{ p,. q,)— 1 
= max{0, 17}— 1 = 16. As the forecast horizon (problem length) increases, the planning horizon increases 
linearly as will be demonstrated by the computational results in the next section. 


6. Computational results 


In this final section we discuss some computational studies of the Forward Algorithm, and compare 
results with those of two specially tailored versions of the simplex method. All computer programs were 
written in FORTRAN, for testing on a time-shared DEC-SYSTEM 2050 at Carnegie-Mellon University. 

We first generated 51 problems with 500 periods, randomly generating demands ( D,), and the initial and 
final work forces. The unit costs were fixed at h = f= 300, u =v = 200 to produce a fixed cost ratio 
K,=(h+ f)/u=K,=(h+f)/v = 3. With these data we incremented the problem length by 25 periods 
starting at 25 and plotted the means of the number of pivots, solution time (msec) and number of horizons 
versus problem length in Fig. 7. Clearly these plots are reasonably linear and by linear regression we find 
for expected values of the three factors: 


xe pivots = 72-825 + 0.8167, Bsstuisontme = ~ 0.852 + 0.2657. 
No horirans = — 0.571 + 0.222T (22) 


where T is the problem length. All the correlation coefficients exceed 0.99. So, for a fixed cost ratio we 
expect the number of pivots, solution time, and number of horizons to be linear in T. The negative 
intercept is due to the start-up of the algorithm. 
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Fig. 7. Means of Number of Pivots, Solution Time (msec), and Number of Horizons versus T, for 51 randomly generated problems, 
solved by the Forward Algorithm. For each problem: (4+ /)/u=(h+ f)/e =3, Wy. W, and demands D, are random. All 
computations performed on a time-shared DEC-SYSTEM 2050. 
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Fig. 8. Mean Solution Time versus T for 15 sets of 21 randomly generated problems solved by the Forward Algorithm. 
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Fig. 9. Mean Number of Pivots for the Forward Algorithm and two versions of the simplex method (see text) versus 7 for 8 sets of 21 
randomly generated problems, one set for each problem length. 
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In Fig. 8 we show a plot of mean solution time (msec) versus T for 15 sets, T= 10. 20, .... 150, of 21 
randomly generated problems. For these data the cost ratios. K, and K,, vary from 2 to 6. Again the results 
are linear in T. 

We next developed a custom tailored simplex code which performed only additions when pivoting. since 
all pivots elements are I’s due to the special structure of the problem. Furthermore, the code does not 
update a row when the entry in the pivot column is 0. We use two entry criteria, the first negative reduced 
cost rule, call this SIMPLEX], and the most negative reduced cost rule, call this SIMPLEX 2. We added 3 
sets of 21 random problems with lengths T = 5, 15, 25 to the data discussed in the previous paragraph, and 
do not consider problems with T > 50. 

In Fig. 9 we plot the mean number of pivots versus T for the Forward Algorithm (F). SIMPLEX 1 (S1) 
and SIMPLEX 2 (S2). All are linear in 7 with correlation coefficients in excess of 0.95. But the mean 
solution times for SIMPLEX 1 and 2 are at least quadratic in T as is shown in Fig. 10. From Figs. 7 and 8 
we observed that the mean solution time for the Forward Algorithm is linear in T. 

From regression analysis we have mean solution times given by 


wy = 0.838 + 0.2337, pe, =0.20277*". py, = 0.0667? *?, (23) 


All have correlation coefficients exceeding 0.98. 
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Fig. 10. Mean Solution Time for two versions of the simplex method versus 7 for 8 sets of randomly generated problems. one set for 
each problem length. 
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Fig. 11. Log—-log plot of Median First Planning Horizon and First Forecast Horizon versus cost ratio K = K,; =(h+ f)/u=K,=(h 
+ f)/v, for 51 randomly generated problems of maximum length 500. 


Finally we investigated the sensitivity of the first planning horizon to the cost ratios K, and K). Initial 
results indicated that the minimum of the two ratios dominates the effect on when the first planning 
horizon occurs. Hence, we set K = K, = K, for these tests. The demand and work force data for the 51 
problems with T = 500 were used here. The cost ratio K was varied logarithmically from 1 to 17.78 after 
which the median first planning horizon was not defined. The results are plotted logarithmically in Fig. 11. 
Clearly the medians for both the first planning and forecast horizons are superlinear in K. Via regression 
analysis we have 


= 0,939K 1-462, = 1,.791K' 488 (24) 


Mrirst Planning Horizon First Forecast Horizon 


with both correlation coefficients approximately 0.95. 

Although Theorem 6 indicates that the Forward Algorithm is no worse than cubic in problem length. 
empirical evidence demonstrates that it is linear while the simplex method is worse than quadratic. In 
addition to its speed, all pivots are performed with comparisons and additions, and memory storage 
requirements are linear in 7. In fact, the solution for time periods through the planning horizon may be 
printed and discarded, freeing computer space for longer problems. The simplex method requires 


J.E. Aronson etal. / A forward algorithm and planning horizon procedure 361 


approximately 1077 storage locations. The Forward Algorithm may be easily executed by hand. augment- 
ing periods as required since no tableau is necessary and degenerate pivots are never executed. 

We conclude that for this specially structured problem, the Forward Algorithm is superior to the 
conventional simplex method. 


Appendix. Proofs of Theorems and Lemmas 


Proof of Lemma 1. For 1 = 1..... 7. set H,= F,=0. W, = W,. If W = D,. then V, = U =0.1f W, < D,, then 
V,=D,- W,. U,=0. If W> D,, then U,= W,- D,. vV,=0. Further, set W,_, = W., and if W. = W, then 
H,.,=F,.,=0, if Wo < W,, then F,., = W,- W.. H,,, = 0. and if W.> W,. then H;,, = W.- W,. 
Fr, =0. 

The above solution is feasible to (P;-), and has at most 2(7 + 1) positive variables. The variables W,, V, 
when W, < D,, and U, otherwise. and H,;, , or F,,, form a basic solution to (P;). Since W,. W.. and D, are 


integer. the solution is integer. 


Proof of Lemma 2. For a degenerate pivot (7) is trivially satisfied. For a nondegenerate pivot, once (7) is 
satisfied by a basic feasible solution, an improved solution can never be attained by violating (7). because it 
would involve firing, hiring and overtime costs, with no benefit. 


Proof of Theorem 1. If (12) is not optimal for (P;-), then there exists an improved solution. But this pivot is 
also possible for solution (11) of (P7 ), which contradicts the optimality of (11). 


Proof of Theorem 2. We abbreviate ‘hiring and/or firing’ as H/F. Suppose we have a basic feasible 
solution to (P,;-). Without loss generality, we consider W,. t= 1,.... T + 1, to be basic. 


“At least two”. 
Case 1, No H/F variable changes. For a non-degenerate pivot, since H, and F, are fixed. by (9), for some 
r, exactly one of V, or U, changes. But by (3), W, must change, a contradiction by (2). 
Case 2. Exactly one H/F variable changes. We consider H, increasing for some ft. Then, to satisfy (2) 
either 
(a) W, increases = W,,, increases = --- = W,.,, increases = (4) is violated, or 
(b) W,_, decreases = W,_, decreases = --- = W, decreases = (2) is violated for r= 1. 
The case for H/, decreasing 1s analogous. The case for F, changing is analogous to that for H, changing. 


“At most two”. 

Let X=({X,...., X,} be the set of as and firing variables in the current basis where X, represents 
either H, or F and i, <i, < Se Sly 

Case 1. Suppose Y,(= 4H, or F,) 1s to enter the basis at the current iteration, and X, leaves the basis. 
Consider i, < ¢. (The case for i, > 1 is analogous.) 

(a) Suppose there exists an X, such that i, <i, <s and X, changes value as a result of the pivot. For 
the pivot to be nondegenerate U.. or V,. but not both, is basic for t’ = 7i,..... 1— 1. Let B be the 
matrix representing the current basis. Since W..... W,.., are basic. we add the columns of B 
corresponding to W,...., Ke to form the echuni vector d. The elements of d are d, =4, = 
Cpa, Apis eS : O otherwise. But @ can be formed by a linear combination of the 
columns in B prec ol s to U, or V,.. t= i,.....1,. X,. and X,. Thus the hypothesized current 
basis is singular, a contradiction. 

(b) Suppose there exists an X, © X, 1 <i,, which is to change value as well as xX, and Y,. Then by 
taking linear combinations of the columns of B as above. we find a linear dependency, a 
contradiction. 

Case 2. Suppose Y,, and H/F variable, enters the basis at a positive value and no X, © X leaves. Then 
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some U,. or V,. must leave the basis. As in Case 1, the columns of B exhibit a linear dependence. 

Case 3. Suppose U,(or V,) enters the basis at a positive value. As in Case 1 the columns of B exhibit a 
linear dependence. 

Because the latter half of the proof depends on existence, the Theorem is true. 


Proof of Theorem 3. It is sufficient to show that A is totally unimodular. This can easily be demonstrated 
by a simple transformation. In A, we add rows T + 2,...,27 + 1 to rowes 2,..., 7 + 1, respectively, to form 
A’. Then by Theorem 3 of [4], A’ is totally unimodular. Therefore, A is totally unimodular. | 


Proof of Lemma 3. We derive condition (i). The others are similar. An additional unit of hiring in period / 
followed by an additional unit of firing in period j affects overtime and undertime only in periods / toj — |. 
k periods will have overtime reduced by one unit, m + n periods will have undertime increased by one unit. 
Thus X(i,/) is the marginal benefit of such a variation in the plan while (A +f) is the extra hiring and 
firing cost to produce this benefit, so (A + f) — X(i, 7) < 0 is the unit cost reduction due to the indicated 
pivot (last column Table 1). Hence the solution is strictly improved by the indicated pivot. 


Proof of Theorem 4. “only if’. Assume S# is an optimal solution to (P;), and that at least one of the 8 
conditions hold. Then, by Lemma 3, a strictly improved solution may be attained by performing the 
indicated pivot, a contradiction. 

“if”. Since (P;) is linear, the simplex method may be employed to find S#. Now, because S} is optimal, 
there exists no nonbasic variable which can enter the basis to produce an improved solution. By Theorem 2, 
a nondegenerate simplex pivot would have to change a pair of H/F variables simultaneously, which is 
possible only if one of the 8 conditions of Lemma 3 is true. 


Proof of Lemma 4. Given overtime and firing in the same period, postpone the firing one period to achieve 
a lower cost leading to a contradiction. Similarly for undertime and hiring. 


Proof of Lemma 5. Either X(r, T+ 1) > X(r+ 1, T+ 1) or else r= T by (15). If D. > W,, apply Lemma 4. 
If D; < W,, then X(T, T + 1) <0. Similarly for the second half of the Lemma. 


Proof of Lemma 6. Follows immediately from the fact that we have augmented an optimal T — | period 
plan. 


Proof of Lemma 7. 
Case (a). 
Dr=W,: X(i,T+1)=X(i,T)-u, 
Y(i, T+ 1)=Y(i, T)-v. 


If any pivoting condition occurs for T, then it should have occurred for T— 1 since X(i, T+ 1) < X(i, T) 
and ¥(i, T+ 1)< Y(i, T), thus contradicting the optimality of S#_,. 
Case (b). 


D,>W,: XG, T+ 1) =X, T) +0, 
Y(i, T+ 1)=Y(i, T)-v. 


In Si, since D; > W;= W,, we have V;= W,—- D,, U-=0, Hr. = Fr., = 0. and W,,, = Wp = W,. 
Since H;.,, = 0 pivots (ii) and (iil) are not possible. 
Since F;,, = 0 pivots (vi) and (vil) are not possible. 


' Special thanks to Egon Balas for this simplified proof. 


J.E. Aronson et al. / A forward algorithm and planning horizon procedure 363 


Since Y(i. T+ 1)< Yi, 7), if either pivoting condition (v) or (viii) occurs, then we contradict the 
optimality of S¥._,. 

Since X(i, T+ 1) > X(i, T), the only possible pivots are (i) and (iv). 

Case (c). 


D,<W. X(i.T+1)=XU.T)—u. 
YU.T+1)=YU. 7) +u. 


The proof is similar to case (b). 


Proof of Theorem 5. (a) Suppose i< q* <i+1, fF >. fF. > 0. and X(i, T+ 1)>0. Now X(i, T+ 1l)= 
X(i, g*) + X(q*. T+ 1), and X(i, q*) < 0 else (iv) would hold for (i, q*). which contradicts the optimality 
of S# ,. Hence X(q*. T+ 1)> 0 as claimed. 

(b) If XG. T+ 1)>(h+/f) then X(r, T+ 1) 2 XU. T+ 1 2(At+f). 

(c) If only pivoting condition (iv) holds, then by (a) it holds for (qg*, 7+ 1). If no V, is driven to 0 as a 
result of the pivot, then X(i. /) remains constant for g* <i <j < T+ 1. hence type (i). (ii), and (iii) pivoting 
conditions cannot occur following the pivot. 

Suppose some Y, is driven to 0 as a result of the pivot. Let ¢, resp. 1, be the index of the minimum resp. 
maximum period for which V, is driven to 0 by the pivot. Since X(/./) decreases by (u +t) for each V, 
driven to 0, i <1 <y by the pivot, then no type (i) pivoting condition is introduced. Since H;, , = 0, pivots 
(ii) and (iii) cannot occur for all (i,/). If pivoting condition (iv) occurred before the pivot for (/, T+ 1) 
with />4,, then this condition remains after the pivot. since X(i, 7+ 1) for i> 1, is unaffected by the 
pivot. 

Next we consider the possible introduction of pivoting conditions (v)—(vili) as a result of the type (iv) 
Pivot. 

If, before the pivot, n(q*. T+ 1, S})=0, that is there were no periods, in which U,= V,=0 for 
q* <t<T, then Y(i,/) remains constant for g* <i <j < T+ 1, hence no type (v)- (viii) pivoting conditions 
are introduced as a result of the pivot. So, we assume n(q*, T+ 1, S;-) = 0 before the pivot. Let 4, resp. t, 
be the index of the minimum resp. maximum period for which V, = U, = 0, g* <1 < T. before the pivot. In 
all four cases, let (i,, /,) be the pair of periods for which the proposed pivoting condition is introduced. 
Note that g* <1, <i,< T+ 1 and i, <2. 

Case 1. Suppose a type (v) pivot is introduced by the type (iv) pivot. If g* </,, then before the pivot 
there existed a type (iv) pivoting condition on(q*. i,) contradicting the optimality of S* |. If g* > i,, then 
before the pivot there existed a type (vi) pivoting condition on (i,. q*), contradicting the optimality of 
STAs 

This and the other three cases are summarized in Table 4. 

(d) Similar to (c). Here some V, is driven to 0. Let 1). f,. ty. 4. 4). 7, be defined as in (c). Note 
r<t,<i,<T+ 1, i, <t,. Since X(i./) either decreases or remains constant for all (/. 7), no type (iv) 
pivoting condition can be introduced. Only for pairs (i, 7+ 1). / > 7, can pivoting condition (i) occur after 
the pivot, else we contradict the optimality of S}. ,. Since there is now hiring in period r we must check 
whether conditions (ii) and/or (iii) are introduced for pairs (i, 7). In fact, this is not possible. Suppose on 
the contrary for some (i, 7), condition (11) now holds, that is X(i. r) > 0. Then this was also true before the 
pivot, so that then X(i. T+ 1)= X(ivr)+ X(r. T+ 1)> X(r. T+ 1). a contradiction. 

Similarly, suppose that condition (ii1) holds. so that F > 0 and X(i.r) > —(A+/). But this was also true 
before the pivot. So X(i, T+ 1) = X(i.r)+ X(r. T+ 1) > 0. Hence condition (iv) held for (i. T+ 1), which 
is a contradiction. As in (c) we summarize in Table 3 the contradictions if we were to assume that a type 
(v) (viii) pivoting condition is introduced, for when n(r. T + 1, S}-) = 0 before the pivot. 

(e) Similar to both (c) and (d). That is to say. no new pivoting condition is introduced if the specified 
pivoting decisions are made. 

(f) By definition w = max{ p q}. The type (iv) pivot does not affect p ( p = 0 still). If g* = q = min{il i € Q) 
and F,. = 0 after the pivot, then q increases, else if F,. + 0 after the pivot. then g is unchanged. 

If no V, is driven to 0 as a result of the pivot. then r is unchanged, since X(/, /) is unaffected by the 
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Table 4 

Summary of the latter part of the proof of Theorem 5(c). Type (iv) pivot 

Proposed qe <a, Existing contradiction 

pivot condition or g* >i, of the optimality of S#_, 

Type (v) q* <4, Type (iv) on (q*.¢,) 
gq’ >i, Type (vi) on (1). q*) 

Type (vi) qe <i, Type (iv) on (q*. 4) 
qe>i, Type (vi) on (14), qg") 

Type (vii) qe <i, Type (iii) on (g*.1,) 
q* >i, Type (vii) on (7). 4%) 

Type (viii) qe <i, Type (iii) on (q*.7,) 
q* >i, Type (vii) on (4). 9°) 


pivot. Suppose some YV, is driven to 0 by the pivot. X(i, /) is unaffected by the pivot for i> 1,, where f, is 
defined as in (c). For i <t, <j, X(i, j) decreases by (u + v) for every period in which V, has been driven to 
0. Hence r cannot decrease. In fact, r either remains constant or increases to that r > ¢, + 1. 

(g) The type (i) pivot sets p = 0. If g > r then it remains constant. If g < r then q is updated so that q>r, 
since q indexes the first firing period following the last hiring period, which is now r. There is some period 
for which V, is driven to 0. Let ¢,, 1, be defined as in (c). X(i,/) is unaffected by the pivot for i> 1t,. For 
i<t, <j, X(i,/) decreases by (u + v) for every period in which V, has been driven to 0. Hence. as in (f),r 
cannot decrease, in fact, r either remains constant or increases so thatr > 1, + 1. 


Proof of Lemma 8. Suppose F,> 0, F..,= ++ = Fuga, =O and s<j, ¥(s.T+1+k)>(h+f); that is, 
suppose we have a type (v) pivoting condition on (s, 7 + 1 +k) where s <j. and of course F, was increased 
as a result of the ‘last’ type (i) or (iv) pivot, before the last k augmentations. Now if a type (v) pivot can be 
performed on pair (s, 7+ 1+), then since Y(s, T+ 1+k)= Y(s,j)+ Yj. T+1+k)>(A+f), by 
definition of s we have Y(s,j)>0, so a type (vi) pivot condition existed on (s./) for S7,,. thus 
contradicting the optimality of ST, ,. 


Table 5 

Summary of the latter part of the proof of Theorem 5(d) 

Proposed r<i, Existing contradiction 

pivot condition orr>i, of the optimality of S$}, 

Type (v) r<i, Type (i) on (r, ¢,) 
r>i Type (v) on (7). 7) 

Type (vi) r<i, Type (1) on (7, 4,) 
r>i, Type (v) on (4). 7) 

Type (vii) r<t Type (ii) on (r.7,) 
r>i, Type (vill) on (4.7) 

Type (viii) r<iy Type (ii) on (r.7,) 


r>i, Type (vill) on (2). 7) 
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Proof of Theorem 6. In the worst case it is possible to perform T type (iv) pivots for each q. or T type (vill) 
pivots for each p (corresponding to X(q. T+ 1)>0. Y(p. 1+ 1)> 0, respectively). This yields at most 
T(T — 1)/2 pivots of these types. Since, for any iteration, if a type (i) pivot is performed. ¢’ <r, and if a 
type (v) pivot is performed, s‘ > s, then for any T period problem, there can be at most T pivots of each 
kind. But for each of these pivots there is at least one less tvpe (iv) or (vill) pivot. Hence the algorithm 
requires at most 7(T— 1)/2 pivots. Since computation on any pivot is O(7). the computation for the 
algorithm is no worse than O(T*). 


Proof of Lemma 9. Since all pivots are defined with respect to p. q. r and s, and since w, r and s for all 
intents and purposes are weakly monotonically increasing. it follows that no change will be made in the 
interval from 1] to 2,. 


Proof of Lemma 10. (a) Suppose to the contrary that r<q,<p< 7+ 1. Now X(r. T+ 1)= X(r.q,) + 
X(q,.p)+ X(p. T+ 1)> X(p, T+ 1), at the same time X(q,. p) < —(A +f), since by optimality, pivoting 
condition (iii) does not hold. Thus X(r. g,) > (A +f). contradicting the fact that pivoting condition (1) does 
not hold. 

(b) Similar to (a). 


Proof of Theorem 7. Obviously w > max{ p,. g,). By definition. at most one of p, and gq, can be positive. A 
type (i) pivot on (7, T + 1) sets p, =r. gq, = 0. p = 0. A type (v) pivot on (s. ¢+ 1) sets gq, = 5. pp, = 0.9 = 9. 
After a type (i) (resp. (v)) pivot, by convention we set s = T+ | (resp. r= T+ 1), since no type (Vv) or (vill) 
pivot is possible for (i, T+ 1+ k) with i< T+ 1. By Theorem 5, r (resp. s) either remains constant or 
increases after the pivot. so that r> max{p,.g,) (resp. s > max{ p,.q,)}). since by definition 2, = 
minfw.r.s)— 1. 


References 


[1] K. Arrow and S. Karlin, Smoothed Production Plans, Studies in the Mathematical Theory of Inventory and Production (Stanford 
University Press, Stanford, CA, 1958). 
{2} M.J. Beckmann, Production smoothing and inventory control, Operations Res (1961) 456-467. 
[3] A. Charnes, J. Dréze and M.H. Miller, Decision and horizon rules for stochastic planning problems: A linear example. 
Econometrica 34 (1966) 307. 330. 
(4] R.S. Garfinkel and G.L. Nemhauser, /nteger Programming (Wiley. New York, 1972) 70. 
(5] F. Hansmann and S. Hess, A Linear Programming approach to production and employment scheduling. Management Tech. 
(January 1960) 46-51. 
[6] A.J. Hoffman and W. Jacobs, Smooth patterns of production, Management Sci. 1(1) (1954) 86 -91. 
(7] C.Holt, F. Modighani, J.F. Muth and H.A. Simon, Planning Production Inventories and Work Force (Prentice-Hall, Englewood 
Cliffs, NJ, 1960). 
[8] M. Klein, On production smoothing, Management Sci. 7 (3) (1961) 286-293. 
(9] H. Kunreuther and T.E. Morton, Planning horizons for production smoothing with deterministic demands: I, Management Sct. 
20 (1973) 110-125. 
{10] H. Kunreuther and T.E. Morton, Planning horizons for production smoothing with deterministic demands: II, Management Sci. 
20 (1974) 1037-1046. 
{11] S. Lippman, A. Rolfe. H. Wagner and J. Yuan, Optimal production scheduling and employment smoothing with deterministic 
demands, Management Sci. 14 (3) (1967) 127-158. 
{12} R.A. Lundin and T.E. Morton, Planning horizons for the dynamic lot size model: Protective procedures and computational 
results, Operations Res. 23 (1975) 711- 734. 
(13) T.E. Morton, Infinite horizon dynamic programming models A planning horizon formulation, Operanons Res. 27 (4) (1979) 
730-754. 
(14) E. Silver, A tutorial on production smoothing and work force balancing. Operations Res. 15 (1967) 985 1110. 


